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Diffuse background of high energy neutrinos arising from interactions of cosmic ray protons with far infrared 
radiation background in extragalactic space is calculated. It is assumed that cosmic ray spectrum at superhigh 
energies has extragalactic origin and is proton dominated.The cosmological evolution of extragalactic sources 
of cosmic ray protons as well as infrared-luminous galaxies is taken into account in the calculation. 

1. Introduction 

The extragalactic background of high energy neutrinos studied in this paper arises from collisions of high en- 
ergy cosmic ray (CR) protons emitted by local sources (e.g., by active galactic nuclei) in extragalactic space, 
with extragalactic background of infrared photons. All previous studies (except the most recent ones 1 1 , 2|) 
took into account the interactions of CR protons with relic photons {T = 2.7 K) only. Evidently, other com- 
ponents of the extragalactic radiation background also must be taken into account in calculation of the extra- 
galactic background of high energy neutrinos. The interval of photon wavelengths which is potentially most 
important from this point of view is (1 — 1000)/im, i.e., the infrared region. At last fifteen years infrared 
astronomy developed very intensively (see, e.g., reviews fSlRl ), and now the rather well-grounded calculation 
of ENB from interactions of cosmic rays with extragalactic infrared photons became possible 1 1 2 1. 

For the concrete calculation of infrared background we chose the paper of Dwek ef a/ Q, where the synthetic 
spectral energy distributions in the large diapason of values of the total luminosities were derived, and the 
simple double power-law form of the local luminosity function (i.e., the present-day number density) of infrared 
galaxies suggested by Soifer et al |6 1. As for the CR spectrum: we do not know well enough the extragalactic 
cosmic ray spectrum and, therefore, we are forced to use the crucial hypothesis about extragalactic origin 
of high energy CR's. Everywhere in our calculations we use the extragalactic model (the crossover energy 
3 X lO^'^eV) and normalize our theoretical CR spectrum on experimental CR data. 

Neutrino flux from interactions of high energy cosmic rays with infrared background was recently calculated 
by Stanev fTl. Our work differs from that of fTl in several respects. The main difference is that in 1 1 1 there 
is no non-trivial cosmological evolution of infrared background: this evolution is assumed to be the same 
as of microwave background radiation. Besides, in Q there is no separating of contributions to ENB from 
infrared and optical diapasons of the background while we calculate ENB from far-infrared part of the radiation 
background only. 



2. Extragalactic CR proton spectrum at different epochs 

To obtain approximate expressions for the spectra and intensities of CR protons in extragalactic space we use 
the cosmological transport (kinetic) equation without integral term, i.e., we work in the continuous energy loss 
approximation introduced, for these problems, by Berezinsky and Grigor'eva Q. This equation is written as 

dn(E,z) d X M iniE.z) 

\ +^ 0^^^ " ' ' ^9{E,z). (1) 
oz obj \ + z 
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Here, n{E, z) is the number density of CR protons with a given redshift z, the function j3{E, z) is the change 
of proton energy in unit interval of z, 

f3{E,z)^j^^^E-tp\E,z)^. (2) 

The first term in rh.s. of eq. (|3 takes into account adiabatic energy losses (those due to the cosmological 
expansion). The function tp^{E, z) is the cooling rate of protons via p7 ttX and P7 pe'^e' reactions at 
the cosmological epoch with redshift z. In our approximation, the cooling rate consists from two parts, 

t;\E, z) = t-^,{E, z) + t-}^f,{E, z), (3) 

(cooling due to interactions with relic and infrared components of the extragalactic radiation background). 

The function g{E, z) in r.h.s. of the kinetic equation describes the combined source of extragalactic cosmic 
rays. This source function can be written in the form 

g{E,z)^p{z)i^{z)f{E)^. (4) 
az 

Here, p{z) is the number density of local CR sources (e.g., AGNs) in the proper (physical) volume, p{z) = 
po(l + is the activity of each local source (the integrated number of produced particles per second), 

ri{z) = (1 + z)"^riQ6{zmax — z). Writing this, we assume that the cosmological evolution of cosmic ray sources 
can be parametrized by power law with the sharp cut-off at some epoch with redshift z,nax (jn and Zmax 
considered as parameters of a model of the combined source). At last, the function f{E) in eq.Q describes a 
form of the differential energy spectrum of the local source. 



3. Far infrared extragalactic background 

For a calculation of the extragalactic radiation background it is convenient to use the cosmological transport 
equation which is analogous to that used in the previous section. The function which must be found is the 
number density of infrared photons at different cosmological epochs, The expression for the 

source function of the kinetic equation is 

g'^'iE^z) = I ^piz,L)S^^{E„L)l- ■ |. (5) 

Here, the function S^^ is the spectral luminosity ("spectral energy distribution") of a luminous local source 
(which is in this case an infrared-luminous galaxy with a given total infrared luminosity L). The function 
p{z, L) in eq.(|5} is the number density of infrared-luminous galaxies with a given luminosity L in the proper 
(physical) volume, which is connected with the local luminosity function p(0, L) by the relation 

p[z,L)=p[Q,j^^^}j{l + zf+-<-, (6) 

where 7^ and 7; are parameters determining the cosmological evolution of the luminosity of the local source, 
and evolution of the comoving density of these sources, respectively. 

The resulting expression for the number density of infrared photons in extragalactic space is (neglecting energy 
losses of these photons during their travelling) is 



3 . / , , ^ ' dt 



n'"iE^,z)= I dx{l±^] [ '^p{x,L)S'^ [E,\±^,l} ^ 



l + x I J L \ 1 + z' J E. 



dx 



(7) 
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Figure 1. a) Proton cooling rates in relic and infrared photon background, b) Extragalactic CR proton spectrum calculated 
with (lower curve) and without (upper curve) taking into account energy losses in infrared background. 



The cooling rate of CR protons in infrared radiation background is expressed through n^^ by the formula 

oo oo 



de- 



(8) 



(7p — E/nip). Here, is the photon energy in the CR proton rest system, (j{er) is the photoabsorbtion cross 
section, /(e, ) is the relative proton energy loss in p7-collision (in the observer system). 



We used the following model for the cosmological evolution of infrared sources: 



Z < Zfiat, 



p{z,L) = p[0 



L 



(1 + ^flat)-" 



{1 + ZflatT''il + Z)^, Zflat<Z<Zc 



(9) 



and p{z, L) — for z > z^ut- The parameters are as follows: 7^; = 3,7; = 2, Zfiat = 1, z^ut — 4. 

On Fig. ^ we show the proton cooling rate on infrared photon gas (together with corresponding function for 
the relic photon case). One can see that at proton energy {3 4) x 10^^ eV, the contribution of infrared 
photons in total cooling rate is noticeable. The results of two calculations of extragalactic CR proton spectrum 
(with and without taking into account energy losses on infrared photons) are shown on Fig. QJ). The theoretical 
CR spectrum at z = was normalized on experimental data; from such normalization we obtained the value 
of the product po?7o entering eq.Q for the source function of the CR proton kinetic equation, « 1.4 x 
lO^^^cTO^'^s^^. For the calculation of the CR spectrum we used the following set of parameters determining 
the spectrum slope and cosmological evolution of CR sources: 7 = 2.5, m = 3.5, z^ax — 5. We assumed that 
extragalactic spectrum of CRs dominates beginning from Eq = 3 x lO^^eV. 

The choice of cosmological evolution parameters jd,Ji, Zfiat, Zcut used in the calculation of infrared back- 
ground and CR spectrum is justified by our prediction of the spectral intensity E^I{E-^) of infrared background 
measured experimentally (Fig.|2t). 
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Figure 2. a) Far infrared radiation background calculated using the model of |5| and eq.j9). b) Extragalactic muon 
neutrino background: "/d = 3, 7; — 2, Zfiat = 1; Zcut ~ 4 (upper curve), Zcut ~ 3 (lower curve). 

4. Extragalactic spectra of high energy neutrinos 

For simplicity we use in this work one-pion approximation, i.e., we assume that in p7-reaction only two 
particles (neutron and pion) are produced. In this case, as is well known the pions have approximately isotropic 
distribution in the center of mass system and, as a consequence, the step-like energy spectra in the observer 
system. The main photoproduction reaction, pc b. + '^infrared 7r+ + n, and subsequent decays of 7r+ and /i+ 
lead to production of v^,i)fj^, v^. In the present work we will be interested only in (j/^ + i>p)-flux. Therefore, we 
add together the neutrino from pion and muon decays. Neutron produced in the photoproduction process will 
decay giving proton, long before the next interaction with photons of the background. Results of the neutrino 
spectra calculation are shown on Fig. |2j). 

5. Conclusions 

In spite of the fact that the density of infrared photons in extragalactic space is much smaller (^1.5 photons /cm^) 
than that of relic microwave photons, the neutrino background appears to be not so small, due to the lower 
threshold for photoproduction and, especially, due to much stronger time evolution of infrared background 
in comparison with that of relic photons. The predicted neutrino fluxes are comparable, more or less, with 
the neutrino fluxes from other extragalactic sources at energy region near 10^^ eV (7-ray bursts, topological 
defects, etc) and deserve further theoretical and, in future, experimental studies. 
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